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The identification of magnetic quantum critical points in heavy fermion metals has provided 
an ideal setting for experimentally studying quantum criticality. Motivated by these experiments, 
considerable theoretical efforts have recently been devoted to reexamine the interplay between Kondo 
screening and magnetic interactions in Kondo lattice systems. A local quantum critical picture 
has emerged, in which magnetic interactions suppress Kondo screening precisely at the magnetic 
quantum critical point (QCP). The Fermi surface undergoes a large reconstruction across the QCP 
and the coherence scale of the Kondo lattice vanishes at the QCP. The dynamical spin susceptibility 
exhibits uj/T scaling and non-trivial exponents describe the temperature and frequency dependence 
of various physical quantities. These properties are to be contrasted with the conventional spin- 
density-wave (SDW) picture, in which the Kondo screening is not suppressed at the QCP and 
the Fermi surface evolves smoothly across the phase transition. In this article we discuss recent 
microscopic studies of Kondo lattices within an extended dynamical mean field theory (EDMFT). 
We summarize the earlier work based on an analytical e-expansion renormalization group method, 
and expand on the more recent numerical results. We also discuss the issues that have been raised 
concerning the magnetic phase diagram. We show that the zero-temperature magnetic transition is 
second order when double counting of the RKKY interactions is avoided in EDMFT. 

PACS numbers: 71.10.Hf, 71.27.+a, 75.20.Hr 



I. INTRODUCTION 

Heavy fermions started out as a fertile ground to study 
strongly correlated Fermi liquids and superconductors^ 
There was a great deal of amazement at seeing a Fermi 
liquid whose quasiparticle mass is over a hundred times 
the bare electron mass; hence the name of the field. It 
was also surprising to find superconductors in an inher- 
ently magnetic environment: the existence of local mag- 
netic moments in these materials is established through 
the observation of a Curie- Weiss susceptibility at inter- 
mediate temperatures (of the order of 100 K), and mag- 
netism is supposed to be "hostile" to superconductiv- 
ity. It was qualitatively understood that the large mass 
reflects the Kondo screening of the magnetic moments, 
which is necessary to overcome magnetism*^ and the 
proper theoretical description of the Fermi liquid state* 
was subsequently achieved. When high temperature su- 
perconductors were discovered in 1986, the development 
of the heavy fermion field was naturally interrupted - at 
least partially. The hiatus proved to be relatively short- 
lived. As it re-emcrged, however, the field acquired a 
considerably different outlook, with the emphasis now 
placed on non-Fermi liquid behavior and magnetic quan- 
tum phase transitions. Since the late 1990s the field has 
become a focal point^ for the general study of quantum 
criticality. The interest in QCPs is by no means unique 
to heavy fermions; it also arises in high temperature su- 
perconductors among other materials^* However, heavy 



fermions are particularly advantageous in one important 
regard. That is, second-order quantum phase transitions 
are explicitly observed in a growing list of this family of 
materials. 

The QCP in heavy fermions typically separates an anti- 
ferromagnetic metallic phase from a paramagnetic metal- 
lic phase. Near the magnetic QCP, transport and ther- 
modynamical properties develop anomalies. The T = 
SDW picture£i2i2ii£ describes the QCP in terms of fluc- 
tuations of the magnetic order parameter - the param- 
agnons - both in space and in time. This theory amounts 
to a <fi 4 theory - with (f> being the paramagnon field - in 
an effective dimensionality of d c s = d + z. Here, z is the 
dynamic exponent, and is equal to 2 in the antifcrromag- 
netic case. So, d e s is above the upper critical dimension, 
4, of the </> 4 theory, for spatial dimensions d > 2. The 
corresponding fixed point is Gaussian. In this picture, 
the non-Fermi liquid properties of the QCP reflect the 
singular scattering of electrons by the paramagnons and 
are unrelated to the process of Kondo screening. 

The most direct indication for the unusual nature of 
the heavy fermion quantum criticality came from inelas- 
tic neutron scattering experiments^*^ The magnetic 
dynamics shows a fractional exponent and uj/T scaling 
over an extended range of momentum space. These fea- 
tures deviate drastically from the expectations of the 
SDW theory. Since the SDW critical fixed point is Gaus- 
sian, spin damping would be controlled by some (dan- 
gerously) irrelevant coupling, and neither fractional cx- 
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ponent nor oj/T scaling would be expected. 

One way to resolve this impasse invokes the break- 
down of Kondo screening at the magnetic quantum crit- 
ical point. On the paramagnetic side, local moments be- 
come entangled with the conduction electrons and, in 
the process, are delocalized and a part of the electron 
fluid. At the magnetic QCP, the magnetic fluctuations 
turn soft and act as a source of dissipation that cou- 
ples to the local moments; this coupling competes with 
the Kondo interactions and destroys the Kondo screen- 
ing. Going from the paramagnetic side to the QCP, the 
electronic excitations depart from those of a Fermi liquid 
and acquire a non-Fermi liquid form. These non-Fermi 
liquid excitations are part of the quantum critical spec- 
trum. Unlike the paramagnons, they are characterized 
by an interacting fixed point. As the control parameter 
is tuned further, into the magnetically ordered side, the 
system is again a Fermi liquid, but the Kondo effect is 
completely destroyed. An important corollary of this pic- 
ture^ is that the Fermi surface has a sharp jump at the 
QCP. 

To microscopically study the magnetic quantum phase 
transition requires an approach that can handle not 
only the heavy fermion and magnetic states but also 
the dynamical competition between the two on an equal 
footing. One suitable approach (if not the only one 
so far) is the extended dynamical mean field theory 
(EDMFT)i 16 i 17 i 18 i 19 

In this article, we discuss the EDMFT studies of the 
Kondo lattice system. In addition to a brief summary 
of the earlier analytical worksji^ we will pay particular 
attention to two issues. The first one deals with the mag- 
netic dynamics near the QCP. A fractional dynamical- 
spin-susceptibility exponent accompanies the destruc- 
tion of Kondo screening^ The second one concerns the 
nature of the zero-temperature transition. When the 
EDMFT is implemented such that there is no double- 
counting of the static component of the RKKY interac- 
tions, the transition turns out to be second orderiSi We 
compare these results with those of some recent related 
worksM^ 



II. THE MODEL, THE PHASES, AND THE 
PHASE TRANSITIONS 

A. Kondo Lattice Model 

We consider a Kondo lattice model, 

H = H f +H c + H K . (1) 
Here, the /-electron component 

u 

describes the interactions between spin-i local moments 
and a = x,y, z are the spin projections. We have taken 



the valence fluctuations to be completely frozen, which 
should be a good description of at least those heavy 
fermion metals that have a heavy effective mass and that 
are undergoing a magnetic quantum phase transition. 
Without loss of generality, we have assumed that a unit 
cell contains one local moment (Si), whose spin-i nature 
reflects the projection onto the lowest Kramers doublet. 
Ifj describes the RKKY exchange interaction between 
the local moments. In the physical systems, the RKKY 
interaction is generated by the Kondo interactions. Here, 
we have taken it as an independent parameter, for two 
reasons. First, it is useful to do so for the purpose of 
specifying the global phase diagram. Second, as it will 
be seen below, in the EDMFT approach to Kondo lattice, 
it is necessary to include this parameter at the Hamilto- 
nian level to treat its effects dynamically (see sections InTl 
and IIVI for details). 

The conduction electron component of Eq. is sim- 
ply 

It is implicitly assumed that the conduction electrons 
alone would form a Fermi liquid, and the residual inter- 
action (Landau) parameters for the conduction-electron 
component alone can be neglected. We will consider the 
number of conduction electrons per unit cell, x, to be in 
the range < x < 1; in addition, we will assume that it 
is not too small, so that the physical RKKY interaction 
between the nearest-neighboring local moments is anti- 
ferromagnetic, and not too close to 1, so that the Kondo 
insulating physics does not come into play. All the phases 
described below are metallic. 

Finally, the local moments interact with s c .i, the spins 
of the conduction electrons, through an antiferromag- 
netic Kondo coupling Jk ■ 

H K = Y, Jk S < - s --<- ( 4 ) 



B. Kondo Effect and Magnetic Quantum Phase 
Transition 

For qualitative considerations, consider the Kondo lat- 
tice model with a fixed value of / and W , with a 
small ratio I/W; here / is the typical (say, nearest- 
neighbor) interaction of Hf, and W the conduction- 
electron bandwidth. We further assume that Hf has an 
Ising anisotropy. In the antifcrromagnctic phase of Hf, 
the spin excitation spectrum is fully gapped. An infinites- 
imal Jk cannot lead to any Kondo screening. Hence, 
the Fermi surface encloses only the conduction electrons, 
whose number is x per paramagnetic unit cell. We label 
this phase AFs- 

On the other hand, when Jk dominates over /, the 
standard Kondo screening does occur. Each local mo- 
ment is converted into a spin-^ charge e Kondo reso- 
nance. The Fermi surface now encloses not only the 
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conduction electrons but also the local moments, the 
total number being 1 + x per unit cell. We label this 
phase PMl- While the existence of this phase is wcll- 
established^ the easiest physical way to sec it is to con- 
sider the limit Jk S> W 3> /. [Since the Kondo state 
restores SU(2) symmetry, we have, without loss of gen- 
erality, taken the Kondo exchange coupling to be spin- 
isotropic in Eq. In this limit, there is a large binding 
energy (of the order — Jk) for a local singlet between S* 
and s ci , and we can safely project to this singlet sub- 
space. In this subspace, x = 1 becomes special: here 
each local moment is locally paired up with a conduction 
electron, and the entire system becomes a Kondo insula- 
tor. For a system of total unit cells, an x < 1 amounts 
to creating (1 — x)Af unpaired local moments, each of 
which is equivalent to creating a hole in the singlet back- 
ground. The Kondo lattice model becomes equivalent to 
an effective single band Hubbard model of (1 — x) holes 
per site, with an infinite on-site repulsion (it is impossible 
to create two holes - there is only one electron in the sin- 
glet to begin with)i2£i2& In the paramagnetic phase, the 
Luttinger theorem then ensures that the Fermi volume 
contains (1 — x) holes or, cquivalently, (1 + x) electrons!. 

These general arguments show that the AFg and PMl 
phases are two stable metallic phases. They differ in 
two important regards. The AFg is magnetically ordered 
while the PMl is not. Equally important, the PMl has 
the Kondo screening while the AFg does not. Increas- 
ing the ratio 8 = Jk/I takes the system from AFg to 
PMl- A key question is this: Does the destruction of 
magnetism and the onset of Kondo screening occur at 
the same stage? If so, the transition is distinctly dif- 
ferent from the T — SDW picture. If not - i.e., if 
the destruction of magnetism happens after the Kondo 
screening has already set in - then the magnetic tran- 
sition can be interpreted as an SDW instability of the 
quasiparticles near the large Fermi surface; the transi- 
tion goes back to the realm of the T = SDW transition 
picture. 

Microscopical studies provide a way to address this 
issue. A suitable method has to capture not only the 
AFs and PMl phases, due respectively to the RKKY ex- 
change interactions and the Kondo interactions, but also 
the dynamical competition between these interactions; 
this dynamical interplay is crucial for the transition re- 
gion. At this stage, the EDMFT method is the only one 
we are aware of which fits this requirement. 

Of course, microscopic studies always have limitations. 
Approximations are inevitably used, in the process of 
solving a Hamiltonian or at the level of the model it- 
self (or both). Controlled approximations, nonetheless, 
provide us with not only ways of understanding experi- 
ments but also intuitions that serve as a basis for more 
macroscopic approaches. This general philosophy is read- 
ily reflected in the EDMFT approach. Even though it 
is "conserving" (i.e., satisfying thermodynamic consis- 
tency requirements), it assumes that the q-dependenecs 
of some irreducible single-electron and collective quan- 



tities \M{uS) and S(w) defined in the next section] are 
unimportant. 

There are several reasons to believe that the EDMFT 
approach is useful for the antiferromagnetic quantum 
transitions at hand. First, because staggered magnetiza- 
tion is not a conserved quantity, the spin damping does 
not have to acquire a strong dependence on q. This is es- 
pecially true in metallic systems, where the dynamic ex- 
ponent z associated with the long-wavelcngth magnetic 
fluctuations is larger than 1. In the SDW case?£ for in- 
stance, z = 2 and M(q, uj) can simply be taken as a 
linear function of \uj\ without any singular dependence 
on q. Whether the q-dependence in M(q, u) is singular 
or regular can be stated in terms of the anomalous spatial 
dimension 77 characterizing the long-wavclcngth fluctua- 
tions in space: a non-zero 77 means that the q dependence 
in M(q, to) is more singular than that [(q— Q) 2 ] already 
incorporated in 7 q . So, if 77 7^ 0, the EDMFT is expected 
to fail, at least for the asymptotic behavior. For instance, 
the classical critical points associated with a finite tem- 
perature magnetic transition in d = 2, 3 must have a 
finite 77; the EDMFT turns out to produce (an artificial) 
first-order phase transition. For a quantum critical point, 
on the other hand, long-wavelength fluctuations occur in 
deft = d+z dimensions. There is then a greater likelihood 
for the vanishing of the spatial anomalous dimension, in 
which case the q-dependence of M(q, u>) is not singular 
and neglecting it will not change the universal behavior. 

Second, as we have already discussed, the different 
classes of magnetic quantum critical points of a Kondo 
lattice can be classified in terms of whether the Fermi 
surface (in the paramagnetic zone), which is large in the 
paramagnetic metallic phase, stays large as the QCP is 
crossed or becomes small by ejecting the local moments. 
Such large vs. small Fermi surfaces are well described 
in terms of whether the Kondo effect is preserved or 
destroyed, which, in turn, are readily captured by the 
EDMFT approach. 

We now turn to the EDMFT studies of the Kondo 
lattice model. 



III. DESTRUCTION OF THE KONDO EFFECT 
WITHIN EDMFT 

A. The EDMFT equations for the paramagnetic 
phase 

The EDMFT approach treats certain intersite (coher- 
ent and incoherent) collective effects on an equal footing 
with the local interaction effects. The EDMFT equations 
have been constructed in terms of a "cavity" method^ 
diagrammatics^i and a functional formalism^ All of 
these constructions yield the same dynamical equations. 
In the diagrammatic language, the EDMFT is seen as 
entirely different from a systematic expansionist in 1 jd 
whose zcro-th order would correspond to the dynamical 
mean field theory (DMFT),^ 2 ^ Instead, the EDMFT is 
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a re-summation scheme that incorporates an infinite se- 
ries of processes associated with the intersite collective ef- 
fects, in addition to the local processes already taken into 
account in the DMFT. Unlike the single-electron proper- 
ties, the collective modes do not have a chemical poten- 
tial. In other words, the bottom of a "band" is important 
and this provides a means for spatial dimensionality to 
come into play in the EDMFT. 

Within the EDMFT, the collective effects are orga- 
nized in terms of an explicit intersite interaction term 
at the Hamiltonian level. For the Kondo lattice model 
described in the previous section, this is the intersite ex- 
change term, Hf of Eq. © . 

There are several ways to see the details of this formal- 
ism. One way is to focus on a spin cumulant^i whose 
inverse, M(w), is colloquially referred to as a spin self- 
energy matrix. While it can be rigorously defined for 
any spinful many-body problem, this quantity is taken 
as q- independent in the EDMFT. 

The dynamical spin susceptibility, on the other hand, 
is q-dependent and is given by 



X a (q,^) 



1 



M a (u;) 



ja 

q 



(5) 



The conduction electron self-energy is still given by 

and the conduction electron Green's function retains the 

standard form, 



G(k,e) 



1 



e + M - £ k - S(e)' 



(6) 



The irreducible quantities, M(uj) and S(e), are deter- 
mined in terms of a self-consistent Bose-Fermi Kondo 
model: 



^pa 



Wimp = Jk S • S e + ^ E p C p(J C p 

,a ~t~ (j)_p a J ~\~ ^ ^ ^p,a tpp^a ^Pp-a * 

(7) 



p. a 



p. a 



The self-consistency reflects the translational invariancc: 



q 

G loc (u) = £]G(k,w) 



(8) 



When combined with the Dyson equations, M a (ui) = 
X » + 1/xg.cM and E(a;) = G^H - l/Gioc(w), 
where Xo.lH = -9 2 Ep 2w Pt a/[uj 2 - w 2 p J and G Q (u) = 
J2 p -Ep), are the Weiss fields, these self-consistency 

equations specify the dispersions, E p and w Pta , and the 
coupling constant g. 

We will focus on the case of two-dimensional mag- 
netic fluctuations, characterized by the RKKY density 



of states p 7 (e) = E q <K e - 7 q ) = (1/27)6(7 - |e|). The 
first of Eq. JHJ becomes 



M» = 7/tanh[7 X ? »] 

= 7+27exp[-27 X ? »] 



(9) 



where the last equality is an expansion in terms of 
exp [— 2Ixi oc (w)], valid when the local susceptibility is 
divergent. 



B. Destruction of the Kondo effect 

Both the Kondo screening and its destruction are en- 
coded in the Bose-Fermi Kondo model, Eq. Q. The an- 
tifcrromagnctic Kondo coupling (Jk) is responsible for 
the formation of a Kondo singlet in the ground state and 
the concomitant generation of a Kondo resonance in the 
excitation spectrum. The coupling of the local moment 
to the dissipative bosonic bath (g) provides a competing 
mechanism. To see this in some detail, we first analyze 
Eq. (J7J alone without worrying about the self-consistency 
conditions. We consider a given spectrum of the bosonic 
bath. 



J~][<5{u - W Pi a) -6({jJ + ' 



J p,aJ 



sgmu 



(10) 



The problem can be studied using an e-expansion^S For 
small g, the Kondo coupling dominates, leading to a 
Kondo screening. A sufficiently large coupling g de- 
stroys the Kondo screening completely, reaching a local- 
moment phase. The transition between these two phases 
is of second-order, and is described by a QCP where the 
Kondo screening is just destroyed and the electronic ex- 
citations have a non-Fermi liquid form. The local spin 
susceptibility has a Pauli form on the Kondo side. The 
destruction of Kondo screening is then manifested in a 
divergent local susceptibility at the QCP. An important 
property that is shared by the Bose-Fermi Kondo model 
with SU(2) spin symmetry, XY spin anisotropy, or Ising 
spin anisotropy, is that Xi oc ( T ) ~ l/ r< S a t the QCP. Here, 
a = x, y, z, a = x, y, and a = z for the SU(2), XY, and 
Ising cases, respectively. 
Correspondingly, 



(11) 



While the critical amplitude, A a (e), depends on the spin 
anisotropy, the critical exponent does not; it is equal to 
1 — e in all cases. 

There is an important point that follows from the 
above analysis which we will use in the following dis- 
cussion of the numerical results. Within the EDMFT 
approach to the Kondo lattice, if the local spin suscep- 
tibility of the Kondo lattice model is divergent at the 
magnetic QCP, the corresponding local problem is sit- 
ting on the critical manifold. In other words, a divergent 
local susceptibility is a signature of the critical Kondo 
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sin 2 (7IX/P) 



FIG. 1: The local spin susceptibility xioc at the quantum 
critical coupling (/ ~ J c ) as a function of the imaginary time 
t. The log of xioc is plotted against the log of sin -1 (ttt//3). 
The long-time limit, r — > /3/2, corresponds to sin _1 (7rr//3) — > 
1. The solid line is a fit in terms of (tt//3)to sin _1 (7rr//3). The 
fit yields a dynamical spin susceptibility exponent a « 0.73. 
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FIG. 2: Plot (log-log) of the local spin susceptibility, xioc vs. 
the Matsubara frequency, u n , at both low frequencies and 
high frequencies. The dotted curve marked "leading order" 
corresponds to a logarithmic dependence of xioc on frequency. 
The dot-dashed curve describes the fitting at high frequencies 
(not shown); the u~ 2 dependence is dictated by the spectral 
sum rule. 



screening and its associated non-Fermi liquid electronic 
excitations. 

These e-expansion results for the Bose-Fermi Kondo 
model were initially used to study the full self-consistent 
problem^ with the self-consistency conditions specified 
in the previous section. The fact that the critical expo- 
nent for the local susceptibility is equal to 1 — e [Eq. 
turns out to be essential for the existence of a local QCP 
solution. The self-consistent solution in the case of two- 
dimensional magnetic fluctuations has e = 1 _ , corre- 
sponding [Eq. ((TIT) ] to a logarithmically divergent local 
susceptibility. 



C. Fractional exponent 

The EDMFT equations in the Ising case (taking only 
the a = z component) were studied numerically in 
Ref. I20H2H using the Quantum Monte Carlo method of 
Grempel and Rozenbergi2i*2£ 

It was found that, at the magnetic QCP, the local spin 
susceptibility is indeed logarithmically divergent. Fig- 
ure ^ shows a log-log plot of the local spin susceptibil- 
ity x\oc{t) vs. sin(7rr//3) at a relatively low tempera- 
ture (T = 0.011T£). It is seen from the figure that 
the zero-temperature limit of the local susceptibility is 
X\oc(t) = A/t. This corresponds to a frequency de- 
pendence that is logarithmically divergent in the low- 
frequency limit. A fit of the data yields the value of the 
amplitude A that is directly related to the critical expo- 



nent of the peak value of the lattice susceptibility. 

Figure |21 shows the logarithmic dependence of xioc(w n ) 
directly, for frequencies smaller than the bare Kondo 
scale Tft. As already mentioned in the previous subsec- 
tion, such a divergent local susceptibility signifies that 
the Bose-Fermi Kondo model is located on the criti- 
cal manifold; correspondingly, there is a destruction of 
Kondo screening at the magnetic QCP of the Kondo lat- 
tice model. 

From the divergent local susceptibility, the self- 
consistency equation 10 also determines the spin self- 
energy and, by extension, the dynamical spin suscepti- 
bility of the Kondo lattice. The inverse peak susceptibil- 
ity, x _1 (Q!^n), where Q is the ordering wavevector, is 
shown in Fig. EI A power-law fit yields a dynamical spin 
susceptibility exponent that is fractional, close to 0.72. 

The fractional critical exponent is only seen at \& n \ < 
Tjf. Likewise, the Fermi- liquid (linear in uj n ) damping 
inside the paramagnetic phase is also seen only at fre- 
quencies up to at most Tj£. 

It is instructive to compare the above results with those 
of Ref. [23, which studied an Anderson lattice model. The 
lowest temperature studied in Ref. |22l is 0.25Tj^ . The 
first non-zero Matsubara frequency, u)\ = 2irT, is already 
larger than T^. As a result, neither the fractional expo- 
nent at the critical coupling I ~ I c nor the linear damp- 
ing at / < I c can be observed. 
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FIG. 3: Inverse peak susceptibility vs. Matsubara frequency, 
at a low temperature (T = O.OllT^-). We have used the 
asymptotic form appearing in the second equality of Eq. JUJ . 
The dashed line is a power-law fitting with the exponent 0.72. 



D. Failure of the local 4 description of the 
Bose-Fermi Kondo model 



quantum critical point is interacting over the entire range 
< e < 1, obeying u/T scaling and hyperscaling. These 
results support the observation^ of u/T scaling in the 
(self-consistent) case of e = 1~. They also imply that 
the Bose-Fermi Kondo model is perhaps the simplest 
model in which the standard description of a QCP - 
in terms of a classical critical point in elevated dimen- 
sions - fails. Physically, the Kondo effect, involving the 
formation of a Kondo singlet, is intrinsically quantum- 
mechanical. In the language of a path integral represen- 
tation for spin, the Kondo singlet formation necessarily 
involves the Berry phase term. It is then natural for the 
destruction of Kondo screening to be inherently quan- 
tum mechanical and, by extension, for the QCP of the 
Bose-Fermi Kondo model to be different from its classical 
counterpart at a higher dimension. A more in-depth un- 
derstanding of the underlying mechanism for this effect 
will surely be illuminating. 



IV. MAGNETIC PHASE DIAGRAM OF KONDO 
LATTICES WITHIN EDMFT 

We now address whether the above results, derived 
from the paramagnetic side, are pre-empted by magnetic 
ordering. To do so, we approach the transition from the 
ordered side. 



It is tempting to consider the Bose-Fermi Kondo model 
as equivalent to a local 4 theory. One maps the Kondo 
coupling to an Ising chain (along the imaginary time r 
axis) with 1/r 2 interactions^ In addition, the bosonic 
bath, with a spectrum of Eq. Q1U[). adds an additional 
retarded interaction, l/r 2_£ . The corresponding local 4 
theory is Z ~ J 2?0exp[— S] where 




(12) 

with the constraint \(f>\ 2 = 1. Indeed, such a local <j> 4 the- 
ory emerges in the large- TV limit of a certain O(N) gener- 
alization of the Bose-Fermi Kondo model. In the N = oo 
case, two of us2S showed that, while the destruction of 
Kondo screening does occur, the fractional exponent is 
absent. Subsequently, Pankov et alw^ demonstrated that 
this conclusion remains valid to order 1/N. 

Recent works have considerably clarified the limita- 
tions of the large N limit and demonstrated the fail- 
ure of the local <fi 4 description of the Bose-Fermi Kondo 
model, at least for e > 1/2 (including the case of the self- 
consistent e = For the local-</> 4 theory, e = 1/2 
is the "upper critical dimension" i2i At e > 1/2, the 
critical point would then be Gaussian, implying viola- 
tions of both lo/T scaling and hyperscaling. A number 
of recent studies on the Bose-Fermi Kondo and closely 
related impurity modcls^SSiSLS, have shown that the 



A. The EDMFT approach from the ordered side 

The EDMFT equations for the antiferromagneti- 
cally ordered phase require normal orderingiSiiLiS of 
Hf- H f = E« hj (| : Sf :: Sf : +{SpSi - ±<Sf)<SJ», 
where the normal-ordered operator is : Sf := Sf — {S"). 
The effective impurity model and the self-consistency 
conditions arc similar to Eqs. and ((HJ, except for the 
following modifications. First, there is a local magnetic 
field - the static Weiss field, h\ oc - coupled to S z . This 
local field, arising through /q, must be self-consistently 
determined by the magnetic order parameter (>5 z )H imp - 

Second, the conduction electron propagators are also 
influenced by magnetism. It turns out that the second 
feature has to be treated with care so that there is no 
double-counting of the RKKY interactions between the 
local moments. In Ref. [^J we avoided double-counting 
the RKKY interactions by working with a featureless 
conduction electron band; in this case, the magnetism 
is driven by the interaction / q already incorporated at 
the Hamiltonian level (in Hf). We will expand on this 
issue in the next Section. 

Our phase diagram is shown in Fig. 21 At I < I c , 
the system is in the paramagnetic metal phase. The co- 
herence scale of the Kondo lattice E* D marks the tem- 
perature/energy below which Kondo resonances are gen- 
erated and the heavy Fermi liquid behavior occurs. In 
particular, the Landau damping is linear in frequency at 
\ui n \ < E* oc . At / > I c , the system has an antifcrromag- 
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FIG. 4: The coherence scale of the paramagnetic heavy- 
fermion phase and the magnetic order parameter, tuaf, vs. 
the tuning parameter, 5 = I/T^. Both quantities are deter- 
mined at T = 0.01 lTjf. Also shown is the Curie constant of 
an unphysical solution without a magnetic order parameter 
at S > S c . Lines are guides to eye. The fact that all three 
curves meet at <5 C implies that the zero-temperature transition 
is continuous. 



netic ground state. There is a finite-temperature first- 
order transition at the Neel temperature, T/y. However, 
Tjv continuously goes to zero as the RKKY interaction is 
reduced. Within the numerical accuracy, it vanishes as 
I — > l£r— the same place where E* oc does so as / — > I~ . 
Plotted in Fig. 0] is the magnetic order parameter, toaf, 
at the lowest studied temperature, T = O.OllT^. Again, 
it is seen to continuously go to zero as / — ► Ij. 

Further support for the second order nature comes 
from the study of a nominally paramagnetic solution at 
I > I c . This solution to the paramagnetic EDMFT equa- 
tions co-exists with the solution to the ordered EDMFT 
equations at I > I c . But this "paramagnetic" solution 
is found to contain a Curie component C /T in the static 
local susceptibility with C = limT->o ^Xioc (T 1 , w n = 0) 
or, equivalently, a jump of magnitude C/T in xi oc {T, u) n ) 
as oj n goes to zero. (Its spin self-energy at zero frequency 
and zero temperature tracks /.) What this implies is that 
the "paramagnetic" solution is not the physical one; in- 
stead, the physical ground state corresponds to the mag- 
netic solution. Nonetheless, the study of this unphysical 
paramagnetic solution is helpful to the determination of 
the zero-temperature transition. It is seen in Fig. |U that 
C extrapolates to zero at the same value of I c (J — ► 1+) 
where E* oc goes to zero (/—►/"). This provides an addi- 
tional consistency check for the critical interaction where 
the paramagnetic phase terminates. 

To summarize Fig. within numerical accuracy, all 
relevant scales vanish simultaneously at I c , and the quan- 
tum transition at zero temperature is continuous. 



B. Avoiding double-counting of the RKKY 
interaction 

To discuss the double-counting issue in some detail, 
we revisit the procedure by which antiferromagnetism 
is treated in the standard DMFTi2& Here, the dynam- 
ical equations are constructed entirely in terms of local 
single-particle quantities; two-particle responses are de- 
rived once the dynamical equations have been solved. On 
the paramagnetic side, the two-particle susceptibility sat- 
isfies the following Bethc-Salpeter equation (in matrix 
form) : 

x _1 (q,^) =Xp\(q,w)-/ ir (q,w). (13) 

Here Xp-h(q, oj) is the particle-hole susceptibility bub- 
ble of the full single-particle propagators G(k, e). The 
triplet particle-hole irreducible vertex has the following 
form (again, in matrix form)j22i2ii 

Jir(q,w)=Xp\ioeM-X^H, (14) 

where Xp-h.ioc(^) is the particle- hole susceptibility bub- 
ble of the full local single-particle propagators Gi oc (e)- 
Combining this with Eq. (|13() implies that, on the para- 
magnetic side, 

X - 1 (q ) cj) = A7 q + x I ocM- (15) 

where 

A/ q = Xp\(q,^-Xp- V o>) (16) 

For our Kondo lattice model, AJ q has the meaning of the 
generated RKKY interaction [after inverting the matrix 
in the (/, c) space]. How can A/ q appear in the (particle- 
hole) spin response and not contribute in the dynamical 
equations for the single-particle quantities? The answer 
lies in the way the Brillouin zone is divided in DMFT 
into "special" q's and "generic" q's4i A/ q of Eq. Ijlbj) 
is non-zero only at "special" q's. The dynamical equa- 
tions are constructed in terms of quantities that are lo- 
cal, i.e. summed over q: the special q's, having mea- 
sure zero, are not important for this summation; only 
"generic" q's have the phase space to contribute to the 
local quantities. To be more specific, consider the hy- 
percubic lattice. AJ q depends on q only through the 

combination Y(q) = (1/d) Yl a =i cos ( l la)^ The disper- 
sion X(q) is 0(1) (in the d — > oo limit) only for "spe- 
cial" q's, e.g. along the diagonals of the Brillouin zone, 
qi = (72 = ■•• = Qd- On the other hand, for "generic" q, 
X(q) vanishes [being formally of order 0(1/Vd), as can 
be seen from the central-limit theorem]. When A"(q) van- 
ishes, Xp-h(q,w) = Xp-h,iocM, so A/ q vanishes! The 
antifcrromagnetic wavevector Q (Q Q = ir for all a), be- 
longs to the set of special q's, so A/q ^ 0. And the anti- 
ferromagnetic instability, from the paramagnetic side, is 
signaled by x _1 (Q,^) = A/ Q + X\ l(u) = 0, at u = 0. 
On the antifcrromagnetic side, the nonzero value for the 
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corresponding A/q, ot , is implemented through the intro- 
duction of the doubling of the conduction electron unit 
cell and different single-particle propagators at the two 
sub-lattices. Formally, this doubling of the conduction 
electron unit cell can be described in terms of an effec- 
tive susceptibility, Xor(Q,w), 

X m -\Q,W) = AJ Q ,or + XtoeM- (17) 

Here, again, A7 Q , or = Xp-h(Q, <*>) - Xp-h,iocH- Thc in ~ 
stability of the ordered state is signaled by Xor _1 (Qj w ) = 
A7q oi . + Xioc( w ) = 0' a ^ so at w = 0. Because the effec- 
tive RKKY interaction incorporated on the ordered side, 
A/q or , is the same as its counterpart on the paramag- 
netic side, A/q, the magnetic transition is in general of 
second order. There is a major limitation to this ap- 
proach. The RKKY interaction, being zero at generic 
wavevectors, does not have enough phase space to dy- 
namically interplay with the Kondo interaction. So the 
self-consistent dynamical equations of DMFT does not 
incorporate A7 q at all, and the Kondo screening is al- 
ways present including at the magnetic QCP. The quan- 
tum critical behavior falls in the SDW category, the same 
as in any static mean-field description of Kondo lattices. 

The EDMFT is introduced precisely to allow this dy- 
namical interplay. Here, an intersitc interaction, as given 
in Hf of Eq. (J2J, is elevated to the Hamiltonian level. The 
Bethe-Salpeter equation (|f 3f) still applies. However, the 
particle-hole irreducible vertex becomes^ 

4r(q,w) = Xp-h,io» - Xfo» - Xo V) - ^ , (18) 

where 7 q is the Fourier transform of the intersitc interac- 
tions, hj , already included at thc Hamiltonian level. We 
have, on the paramagnetic side, 

X -1 (q,w) = A/ q + / q + MH , (19) 

where M(uj) = Xioc( w ) + Xo 1 ( LJ ) is the spin self-energy. 
Likewise, we can write the effective susceptibility that 
comes into the stability analysis of the ordered phase as 

Xor _1 (q,w) = A7 q , or + 7 q + M(u) . (20) 

It was shown in Ref. that the EDMFT can be rig- 
orously formulated only when all q are considered to be 
generic. [Otherwise, 7 q is formally of order 0{d) at the 
special q's, and no paramagnetic phase could exist.] This 
implies that A7 q = for all q. From the Kondo lattice 
point of view, this is equivalent to saying that we will only 
use 7 q to represent the RKKY interaction and will not in- 
corporate additional, generated RKKY interactions from 
the fermion bubbles (illustrated in Fig. 7 of Ref. 0). 

In order to be consistent, one would also need to de- 
mand that A7 q or = on the ordered side. Otherwise, we 
would be counting additional contributions to thc RKKY 
interaction on the ordered side that were absent on the 
paramagnetic side. This requirement (A7 qor = 0) was 
achieved in Ref. |^] by using a featureless conduction- 
electron band. The latter ensures that all wavevectors 



are generic in thc sense defined earlier. Within this pro- 
cedure, the magnetic ordering is entirely driven by 7 q , 
and the instability criteria from the paramagnetic and 
ordered sides coincide. Therefore, the quantum tran- 
sition is naturally of second order, as was indeed seen 
numerically in Ref. |2ll see Fig. 0] above. 

The procedure used in Refs. and amounts, in 
our language, to keeping A7q = while A7q jOT ^ 0. 
On the paramagnetic side, all wavevectors, including the 
ordering wavevector Q, are taken as generic, and A7q = 
0, as in all EDMFT schemes. On the ordered side, Q is 
considered as one of the special wavevectors in the sense 
defined earlier and, through thc conduction electron unit- 
cell doubling, A7 Q , or ^ (as in DMFT). The ordered 
side then has an added energy gain, and the magnetic 
quantum transition is of first order. (That an EDMFT 
approach to Kondo lattices which incorporates a DMFT- 
likc fermion bubble on the ordered side alone yields a first 
order transition at zero temperature was in fact already 
recognized in2i.) The procedure would actually lead to 
a first-order magnetic transition at zero-temperature in 
any itinerant model, including any T = SDW transition 
without any Kondo physics. 

We close by noting that the different EDMFT schemes 
that we have discussed can be equivalently seen as differ- 
ent local approximations to a Baym-Kadanoff-type func- 
tional. 



V. EXPERIMENTS AND OTHER 
THEORETICAL APPROACHES 

An important manifestation of the destruction of the 
Kondo screening is that /-electrons participate in the 
Fermi volume on the paramagnetic side but fails to do 
so on the antiferromagnetic side. There is a sudden re- 
construction of the Fermi surface across the magnetic 
QCP. Fairly direct electronic evidence for this effect 
has appeared in thc recent Hall-cffcct measurement in 
YbRli2Si2i^ The Hall coefficient shows a rapid crossover 
as a function of the control parameter — magnetic field 
in this case. The crossover sharpens as temperature is 
lowered, extrapolating to a jump in the zero-temperature 
limit. The jump occurs at the extrapolated location of 
the magnetic phase boundary at zero temperature. Re- 
lated features have also been observed in YbAgGe4£ 

A more direct probe of the Fermi surface comes from 
the de Haas- van Alphcn effect. Recent dHvA mcasure- 
mem™ in CcRhIii5 provides tantalizing evidence for a 
large reconstruction of thc Fermi surface, with a diver- 
gent effective mass, at a QCP. Specific-heat measurement 
under magnetic field^ points towards the possibility that 
CeRhlns undergoes a second-order magnetic quantum 
transition at the magnetic field of the strength used in 
the dHvA experiment. If the existence of the magnetic 
QCP is indeed established, CeRhlns will provide more in- 
sights into quantum criticality than CeRli2Si2. In thc lat- 
ter system, a large Fermi-surface reconstruction has also 
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been seen in the dHvA measurements ^ but the zero- 
temperature transition is likely to be first order with a 
large jump in the magnetic order parameter across the 
transition. 

The fractional exponent and uj/T scaling in the mag- 
netic dynamics have been seen, since early on, in the anti- 
ferromagnetic QCP of Au-doped CeCugAA^ (whose mag- 
netic fluctuations have a reduced dimensionality) and in 
some frustrated compounds (On the other hand, the 
SDW behavior is observed in the magnetic dynamics of 
Ce(Rui_ x Rh x )2Si2r^ which has quasi-3D magnetic fluc- 
tuations.) Related non-trivial scaling exponents - that 
are relatively easy to connect with theory - have come 
from the Griineisen ratio 4& 

Theoretically, there have also been efforts to study 
the Kondo lattice systems using certain mixed-boson- 
fermion representations for the local-moment spin opcra- 
tors42i£&i&L Such auxiliary-particle representations set up 
the basis for a picture with spin-charge separation. How- 
ever, it has been hard to use this formalism to properly 
capture the Kondo-scrccncd Fermi liquid phased mak- 
ing it difficult to study its destruction as well. 

It may also be possible to describe the destruction of 
Kondo screening in terms of the static mean field the- 
ories based on slave boson and an RVB order parame- 
ter, supplemented by gauge-field fluctuations. The cor- 
responding phase diagram has recently been studied in 
some detail^ The magnetic transition and destruction of 
Kondo screening are found to occur at different places in 
the zero-temperature phase diagram^ so the magnetic 
quantum transition is still of the SDW type. We believe 
that this is a reflection of the static nature of the mean 
field theory. 

Finally, it is instructive to put in the present context 
the QCP proposed for the transition from an antiferro- 
magnet to a valence-bond solid in frustrated quantum 
magnets^ Dubbed a "deconfined" QCP, it has certain 
properties that may be qualitatively compared with the 
local quantum criticality: the QCP - containing exotic 
excitations - is surrounded by two conventional phases, 
and the corresponding energy scales of both vanish as the 
QCP is approached. Hence, it would be enlightening to 
explore the concrete connections (if any) of this approach 
with the physics of the destruction of Kondo screening. 
For this purpose, it would be necessary to cither con- 
struct microscopic spin models for the deconfined QCP 
or reformulate the Kondo screening beyond microscopic 
approaches. 



VI. SUMMARY AND OUTLOOK 



To summarize, we have discussed some of the micro- 
scopic approaches underlying the local quantum crit- 
ical picture. Beyond the initial studies based on an 
e— expansion renormalization group method, the most 
extensive investigations have been carried out in Kondo 
lattice models with Ising anisotropy. The latter have al- 
lowed the study of both the destruction of Kondo screen- 
ing and the concomitant fractional exponent and lj/T 
scaling in the magnetic dynamics. We have also discussed 
the magnetic phase diagram and summarized the evi- 
dence for the second order nature of the magnetic quan- 
tum phase transition. The EDMFT studies of Kondo 
lattice models with continuous spin symmetry [SU (2) or 
XY] are mostly confined to the e-expansion studies. Ef- 
forts to access the quantum critical point in these sys- 
tems, beyond the e-expansion, are still underway. A dy- 
namical large-A^ limit, for instance, has recently been 
shown to be promising^ 

The microscopic approaches described here have shown 
that critical modes beyond the order parameter fluctu- 
ations exist, and the modes are associated with the de- 
struction of Kondo screening. These insights have a num- 
ber of phenomcnological consequences - not only in mag- 
netic dynamics but also in the Fermi surface properties 
and in thermodynamics - which have been supported by 
the existing and emerging experiments. The insights will 
also help the search for the field theory that describes 
quantum critical heavy fermions. Finally, they may very 
well be broadly relevant to the exotic quantum critical be- 
havior in other strongly correlated systems such as doped 
Mott insulators. 
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